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ABSTRACT 

We investigate the collapse and fragmentation of low-mass, trans-sonically turbulent 
prestellar cores, using SPH simulations. The initial conditions are slightly supercritical 
Bonnor-Ebert spheres, all with the same density profile, the same mass {M^ = 6.1 M©) 
and the same radius {R^ = 17, 000 AU), but having different initial turbulent velocity 
fields. Four hundred turbulent velocity fields have been generated, all scaled so that 
the mean Mach number is = 1. Then a subset of these, having a range of net 
angular momenta, j, has been evolved. The evolution of these turbulent cores is not 
strongly correlated with j. Instead it is moderated by the formation of filamentary 
structures due to converging turbulent flows. A high fraction 82%) of the protostars 
forming from turbulent cores are attended by protostellar accretion discs, but only a 
very small fraction {^-^ 16%) of these discs is sufficiently cool and extended to develop 
non- linear gravitational instabilities and fragment. 
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1 INTRODUCTION 

One of the key unsolved problems in star formation is 
how a low-mass prestellar core is converted into a small 
stellar system, and what role is played in this pro- 
cess by protostellar accretion discs. This problem is cru- 
cial to understan d ing how the core mass function (e.g. 
Motte et ai] 1 19981: iTesti fc Sargmtl 1 19981: I Johnstone et al.l 



sembly of the molecular cloud is dissipated as it cascades 
to smaller scale^. Thus the larger, more diffuse structures 
within molecular clouds, those that form star clusters, have 
supersonic t urbulence. Thi s is th e regime t hat ha s been 
simulated by iBonnell et all l|2008i ) and lBatei (\2009ih . Con- 
versely, the smaller, denser prestellar cores, which form sin- 
gle stars and small stellar systems, display trans-sonic or 



200ll : iNutter fc Ward-ThompsonI l2007l : lAlves et al.l |2007| ) 
maps into the system and ste llar initial mass functions 
Kroupa| '2001'; ' Chabriej |2003i ). In a companion paper 
Walch et al. 200cr hereafter W09) we have simulated the 



ispl: ., 

subsonic turbulence (e.g. Myers et al.l 0*9911 : lAndre et al.l 



evolution of rigidly rotating cores having different net an- 
gular momenta. Such cores collapse to produce a primary 
protostar surrounded by a protostellar disc, but only cores 
with high angular momenta produce protostellar discs which 
are sufficiently cool and extended to fragment and thereby 
spawn secondary companions. In this paper we simulate 
the collapse and fragmentation of mildly turbulent prestel- 
lar cores having the same range of net angular momenta as 
those simulated in W09. These mildly turbulent cores tend 
to fragment in a completely different way: first a filament 
forms, and then protostars precipitate out of the filament. 

The context for this study is the paradigm in which 
molecular clouds are assembled, form stars and dis- 
perse, on a dynamical timescale (e.g. ElmegreerJ 2000l : 

' mi 



1 19931 : Ib arranco fc Goodmanlll998l ). This is the regime with 
which the present paper is concerned. 

Previous numerical work on the collapse and fragmenta- 
tion of turbulent low-mass cores has concentrated on the sta- 
tistical properties of the resulting protostars (i.e. masses and 
multiplicities) rather than the phenomenology of fragmen - 
tation. For example, iGoodwin et~all (|2004bl : l2004d : |2006| ) 
follow the collapse and fragmentation of 5.4 M© cores hay- 
ing Plummer-like density profiles and very low initial levels 
of turbulence characterised by 7TURB =0.05, 0.10, and 0.25, 
where 



Tturb 



(1) 



IVazguez-Semadeni et al.]|2000l : iBurkert fc Hartman nf l2004l ) 
In this paradigm the turbulent energy delivered by the as- 
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tional energies. Goodwin et al. use SPH with sink particles 

^ Additional turbulent energy may be injected by radiative and 
mechanical feedback from stars, but Larson's scaling relations 
(as discussed in Section 2.2) strongly suggest that most turbulent 
energy is injected on the largest scales and cascades to smaller 
scales 
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and a barotropic equation of state, and perform multiple 
realisations for each value of 7turb • average these sim- 
ulations produce one star per initial Jeans mass, and 80% 
of cor es spawn multiple systems. Recently, lAttwood et alj 
have repeated these simulations, solving the en- 
ergy equation explicitly, and using an approximate method 
to treat the associated transfer of cooling radiation (see 
IStamatellos et a"l]|2007f ). These simulations tend to produce 
more protostars - and in particular more very low-mass 
protostars - than the earlier simulations performed using 
a barot ropic equa tion of state. 

O ffner et al.l (2008) model the collapse of turbulent 
cores, using AMR and a barotropic equation of state. Their 
cores are produced in a large-scale simulation as part of a 
molecular cloud, with either driven or decaying turbulence. 
Individual cores with masses and sizes comparable to those 
treated here are then followed at higher resolution; this has 
the advantage that the cores have consistent velocity and 
density fields from the outset. Their principal finding is that 
simulations with decaying turbulence form more low-mass 
protostars than simulations with driven turbulence, because 
the decline in turbulent support allows forming protostars 
to fall into the centre of the core where they then experience 
competitive accretion and dynamical ejections. 

The plan of this paper is as follows. In Section 2 we 
describe the initial conditions, the numerical code and the 
constitutive physics; we concentrate on the motivation for, 
and the generation of, the initial velocity field, since this is 
the only aspect that is different from W09; more detail on 
the other aspects can be found in W09. In Sections 3 through 
5, we present the results, and discuss the role played by fila- 
ments and discs in the formation and growth of protostars. 
In Section 6 we summarise our main conclusions. 



2 INITIAL CONDITIONS, NUMERICAL 

METHOD, AND CONSTITUTIVE PHYSICS 

2.1 The initial core density profile 

As in W09, we model cores as Bonnor-Ebert spheres trun- 
cated at dimensionless radius — 6.9, and then we increase 
the density everywhere by 10%. The gas is pure H2 at 28 K, 
so the isothermal sound speed is Oq = 0.34kms"^ The 
central density is = 10~^*gcm~^, so a core has mass 

~ 17,000AU, boundary 
^ K, and ratio of thermal 



M, 



e.iMp 



Q _ ^^.j_i.j.0, initial radius 7?^ 
pressure -P^xt — 5-5 x 10^ cm 
to gravitational energy a — 0.74. The freefall time at the 
centre is 67 kyr and the sound crossing time is 250 kyr. 



2.2 The initial velocity field 

The velocity fields in molecular clouds appear to subscribe 
(with some scatter) to a scaling law of the form (j(A) oc A"*, 
where A is the linear size of a coherent feature, and (t(A) 
is its internal velocity dispersion. If these velocities derive 
fro m turbulen ce with a powe r law Pk oc fc" , then n — —3 — 
2q (|Mvers fc Gammic 1999). iLarsoiil l|l98lf ) finds q ~ 0.38, 
corresponding to n = —3.76 (close to the n = 11/3 = 3.67 
of Kolmogorov 1936). iGoodman et all (|1998D find q ~ 0.5 
on scales A > 0.1 pc, corresponding to n = —4. 

Cores created in a turbulent molecular cloud 
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If 


1.03 


0.393 
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148 


0.38 


0.42 


213 


2b 


1.37 


0.466 


137 


141 


0.40 


0.70 


1572 


3a 


1.71 


0.327 


126 


136 


0.33 


0.62 


274 


4b 


2.59 


0.630 


159 


180 


0.26 


0.40 


250 


6d 


2.72 


0.483 


112 


138 


0.22 


0.85 


676 



Table 1. Column 1 gives the ID of the simulation; the number 
corresponds to the ID of the rigidly rotating simulation having the 
same net angular momentum, and the lower-case letter identifies 
which realisation of the velocity field has been adopted. Column 
2 gives the mean specific angular momentum, j. Column 3 gives 
the ratio of rotational to gravitational energy, 7turb . Columns 4 
and 5 give the times at which the primary protostar forms and the 
simulation is terminated. Columns 6 and 7 give the final masses 
of the primary protostar and attendant disc. Column 7 gives the 
radial extent of the disc. 



will have irregular internal velocity fields, and in 
general they will in he rit a n et angular mo m entum 
JColdsmith fc ArguiUal 1 19851: 
iGoodman et al.l Il993[: | Barranco fc GoodmanI 



CascUi ct all 120021 ). Ijiiina et all l| 19991 ') 



1995 



1998 



estimate that 

a low-mass core typically has a ratio of turbulent to 
gravitational energy in the range < Itudb ~ 0-5, 
and mean specific angu la r mo mentum j ~ lO^^cm^s^^. 



iBurkert fc Bodenheimeil l|2000l ) have shown that these 
features can be reproduced if the turbulence has n = —3 or 
n = -4. 

We have therefo re created turbulent veloc i ty fiel ds with 
the same ansatz as iBurkert fc Bodenheimerl (|20oO). First 
we generate a random Gaussian velocity field with oc 
k~*, populating the wavenumbers kj^^j^ ^ fc sC ^maxi with 
'^MiN = 1 (corresponding to the diameter of the cloud) or 
fc^^jjj — 2 (corresponding to the radius of the cloud), and 4 ^ 
'^MAx ^ 8. Then we map this field onto a uniform 128^^ grid, 
covering the whole computational domain. Next we scale the 
velocities of the grid points so that their mean Mach number 
is unity (i.e. the turbulence is transonic). Finally, we obtain 
the velocities of individual SPH particles by interpolating 
on this grid. 

Four hundred velocity fields have been generated in this 
way. We have then selected a subset of these velocity fields 
having the same - or very nearly the same - net angular 
momenta as the rotating cores simulated in W09. Key pa- 
rameters of this subset are presented in Table 1; these are 
the cores whose evolution we have simulated, and whose 
evolution we discuss in this paper. We find that only when 



we use k. 



1 (i.e. there is turbulent energy on the scale 



of the diameter of the core), does the resulting distribu- 
tion of core angular momenta match the observations, with 
j ^ lO^^cm^s"^. If we set fc^j^ = 2 (i.e. the largest tur- 
bulent wavelength is the core radius), then the mean spe- 
cific angular momentum is an order of magnitude lower, 
j ~ 10^" cms-^ 
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2.3 The VINE code 



3 OVERVIEW OF SIMULATIONS 



As in W09, the simulations have been performed with 
the Tree-SPH c ode V INE VINE is describe d fully in 
IWetzstein eFal] (|2008D and iNelson et al.l (|2008D . It is par- 
allelised with OpenMP directives. It invokes a leapfrog 
integrator, and individual particle time steps. We adopt 
a CFL tolerance parameter of 0.1. Gravitational accel- 
erations ar e estimated using a tree, with opening angle 
e = 0.005 (|Springel et all 1200 ll ). The gravitational soften- 
ing length is set to th e hydrodynamical smoothing length 
l|Bate fc Burkertlll997l ). which is adapted so that each parti- 
cle has A/'i^EiB = 50 ± 20 neighbours. Hydrodynamical forces 
are treated with periodic boundary conditions, but gravita- 
tional forces are not. Art ificial viscosity is treat e d usin g the 
standard prescription of iGingold fc MonaghanI l| 19831) with 



1 and f3^ 



2, plus the Balsara switch ( Balsaral 



1 19951 ) to ensure a better treatment of shear flows. 



2.4 Equation of state and molecular line cooling 

As in W09, we compute the gross thermodynamics on the as- 
sumption that the gas is pure , with ratio of specific heats 
7 = 7/5 and molecular weight /i = 2. We solve the energy 
equation, with the following prescription for the radiative 



coohng. At low densities, p < 



10" 



gem ^, we use 



the cooling rates computed bv iNeufeld et all |l993), with 
the constraint that the temperature is not allowed to fall 
below T = 9K. At high densities, p > PcpuT, we switch off 
the radiative cooling, and the gas then evolves adiabatically 
(b ut not isentropically) . The same procedure was adopted 
bv lBaneriee et all l|2004 ). 



2.5 Resolution 

As in W09, the core is modeled with 430,000 SPH parti- 
cles, each having mass m^p^ = 1.4 x 10~^Mq. Hence the 
minimum resolvable mass is 2Afj^^j^^mgp^ = 1.4 x 10"'' Mq 
and the Jeans mass is always resolved. In addition there are 
24,000 ambient particles exerting the external pressure that 
contains the core. We do not use sink particles. Instead we 
impose a minimum smoothing length, /i^m = 2 AU. 



2.6 Definitions 

As in W09, the protostars in these simulations comprise all 
material having density p > 10~^^ gem"'', and the first pro- 
tostar to form is the primary protostar. Similarly the atten- 
dant protostellar discs comprise all material h aving density 
in the range 10"^^ gcm"^ < P < 10~" gcm~^: IOffner et all 
l|2009l ) use a similar definition for discs. Filaments are then 
defined as all material with density p < 10~^^gcm~^ and 
temperature T > 20K. The material in a filament is warm, 
because it is heated in an accretion shock as it enters the 

X + y + z ) and 



Lo = (^x^ + j/^)^^'^ to distinguish, respectively, distance from 
the centre of co-ordinates and distance from the z-axis. 
When analysing protostellar discs, we shift the co-ordinate 
system so that the protostar is at (x, y, z) = (0, 0, 0) and z 
points along the disc's shortest axis of inertia. 



THE EFFECT OF THE TURBULENT VELOCITY FIELD. In turbu- 
lent simulations, regions of large | V ■ v| are present through- 
out the volume of the core, from the outset, and the spherical 
symmetry of the core is therefore immediately broken. This 
means that at the start of the simulation, some parts of the 
core are expanding (V ■ v > 0) and other parts are con- 
tracting (V ■ V < 0). Some of the expanding parts may even 
expand forever and disperse, but most remain bound and 
eventually fall back towards the centre of mass. The con- 
tracting parts tend to form coherent filamentary structures, 
and once these filaments become self-gravitating, protostel- 
lar fragments condense out of them. Sometimes two neigh- 
bouring condensations form from the same filament and sub- 
sequently merge to form a single protostar. On other occa- 
sions two condensations form at well separated locations, 
and in this case they are likely to condense separately and 
form a wide binary protostar. However, such wide binary 
protostars appear to be quite rare. 

THE LOCATION OF PROTOSTAR FORMATION. Protostel- 
lar condensations are not in general centrally located within 
their birth core, and in one simulation a protostar actu- 
ally forms outside the initial boundary of the birth core 
(Run 4b). Typically a protostar has a peculiar velocity of 
~ O.lkms"^, relati ve to the birth c ore. T his is consistent 
with the results of iBeichman et al] (|l986l ). who looked at 
the locations of IRAS sources in the vicinity of dense NH3 
cores having comparable masses and radii to those simulated 
here. About 40% of the IRAS sources had visible counter- 
parts (e.g. T Tauri stars) that lay just outside their associ- 
ated core, implying either that they formed there, or - more 
probably - that they migrated there. During a mean lifetime 
of lO'^ yr and given a mean peculiar velocity of O.lkms"^, 
most mature T Tauri stars would be found just outside their 
birth cores. 

PROTOSTELLAR DISCS. A condensation is continually 
fed material from the filament in which it is embedded. 
Within the condensation, the low angular momentum mate- 
rial rapidly contracts to form a protostar, and the remaining 
material orders itself into a protostellar disc and then slowly 
accretes onto the protostar, on a timescale determined by 
the strength of the gravitational torques which redistribute 
angular momentum within this disc. These protostellar ac- 
cretion discs tend to be relatively compact and hot, and 
hence stable against fragmentation. Because accretion onto 
the discs is lumpy and irregular, they are continually dis- 
turbed. Consequently the discs do not settle down into a thin 
quasistatic equilibrium and there is strong vertical mixing. 

TUMBLING FILAMENTS. The angular momentum of an 
accretion disc is determined by the dynamics of the filament 
out of which it condenses. If the filament is tumbling (i.e. 
spinning about an axis at a sigiuficant angle to its length), 
the material flowing into the protostellar condensation at 
later times tends to have higher speciflc angular momentum, 
and as a consequence the disc will become quite extended, 
and may fragment. However, this is a rare occurrence, be- 
cause the fllaments do not often tumble signiflcantly; their 
orientation in space tends to be approximately constant. 

TWO REPRESENTATIVE EXAMPLES. In order to illustrate 
these properties, we describe in detail two runs from the en- 
semble of simulations we have performed. Runs lb and 6d. 
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logio S [g/em^] 



logio S [g/em^] 




-6000 6000 



Figure 1. False-colour column-density images for Runs lb and 6d, with column-averaged velocity vectors super- 
imposed. The three lefthand columns show frames from Run lb; successive rows correspond to the times t = 
3.94, 15.76, 39.39, 78.79, 118.2 and 152.3 kyr. The three righthand columns show frames from Run 6d; successive rows correspond 
to the times t = 1.97, 9.85, 19.70, 39.39, 99.07 and 140.8 kyr. The density scale is shown at the top of the page, and a lkms~^ velocity 
vector is given on the a; = slices. All frames are (24000 AU)^. 



Protostellar discs formed from turbulent cores 5 



We divide the detailed discussion into two section. Section 4 
deals with the dynamics of filament formation and the for- 
mation of condensations. Section 5 deals with the dynamics 
inside a condensation, and in particular the properties of 
protostellar accretion discs. 

THE INFLUENCE OF THE NET ANGULAR MOMENTUM. Al- 
though these simulations represent the extremes of the dis- 
tribution of specific angular momenta, j (see Table 1), we 
stress that j is not the only - nor even the main - factor 
that leads them to follow distinctly different evolutionary 
paths. Nor is this due to different levels of turbulence, since 



the two simulations have similar values of 7^ 



The diver- 



gence is mainly attributable to the stochastic details of the 
initial turbulent velocity fields. In Run lb, the initial velocity 
field is mainly expansive, and so the primary protostar only 
forms after this expansion has been reversed and the mate- 
rial falls back on itself. The accretion disc also grows very 
slowly, so it is never sufficiently massive and/or extended to 
fragment. Conversely, in Run 6d the initial velocity field is 
compressive, and the primary protostar forms more quickly, 
although its condensation is then held up somewhat by rota- 
tion, and the subsequent growth of the primary protostar by 
accretion from the attendant protostellar disc is slow (much 
slower than in Run lb). 



4 GLOBAL DYNAMICS AND FILAMENT 
FORMATION 

The formation of filaments. Figure [T] shows false-colour 
column-density images on the principal cartesian planes for 
Runs lb and 6d, with column-averaged velocity vectors su- 
perimposed. In Run lb, (the three lefthand columns of Fig. 
[T]), the initial velocity field has large positive V ■ v and this 
breaks the cloud into two pieces. When these two pieces turn 
around and fall back together, they form a filament, and the 
primary protostar condenses out towards one end of this fil- 
ament. Since the filament is not tumbling significantly, the 
flux of material approaching the primary protostar carries 
little angular momentum, and it forms a rather compact 
hot disc, which is unable to fragment. By contrast, in Run 
6d (the three righthand columns of Fig. [!}, the initial ve- 
locity field is both compressive and rotational. As a result 
it quickly produces a tumbling filament, and then a proto- 
star condenses out near the centre of the filament. Because 
the filament is tumbling, the fluxes of material approach- 
ing the protostar from opposite ends of the filament become 
offset from one another, and therefore deliver high angular 
momentum material into the condensation; consequently its 
disc is quite extended and prone to gravitational instabili- 
ties. 

CORE DENSITY PROFILES. Fig. [2]shows p(r), the spheri- 
cally averaged density profile, at various times from near the 
start of the simulation to just after the primary protostar 
forms at t^. In Run lb, the primary protostar condenses 
out near r = 15, 000 AU (quite close to the edge of the ini- 
tial Bonnor-Ebert sphere). In 6d, two condensations form 
out of the same filament, one at r = 2, 000 AU that becomes 
the primary protostar, and one at r = 11, 000 AU; the latter 
may eventually become a wide companion, but the simula- 
tion had to be terminated before it had condensed out. 
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Figure 2. The mean (spherically averaged) den- 
sity at selected times, (a) Run lb at t — = 
-152, -126, -100, -75, -49, -24 and kyr. (b) Run 6d 

at t - to = -102, -83, -63, -43, -24, -4 and + 5 kyr 
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Figure 3. Distances from the centre of the core at which primary 
protostars form vs. initial turbulent energy content 7rrurtB °f ^he 
core. The shaded bar marks the transition region defined in the 
text, and the dashed line is the edge of the initial Bonnor-Ebert 
sphere. 
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Figure 6. The mass of the primary protostar as a function of 
time, for all the turbulent simulations. 



THE TRANSITION ZONE. In the central zone of a Bonnor- 
Ebert sphere — dln[p]/dln[r-] ~ 0, and in the outer zone 
— dln[p]/dln[r] ~ 2. We therefore define a transition zone 
where 0.5 ^ — dln[p]/dln[r] ^ 1.5. In terms of the di- 
mens ionless isothermal function ijChandrasekhar fc Ware3 
1 19491 ). the transition zone corresponds to 0.5 ^5 ^i/)' ^ 1.5, 
and hence 1.32 ^ ^ ^5 2.90. In our simulations, this trans- 
lates into 3,250AU r 7, 150AU. Fig. [3] shows 7turb 
plotted against the distance from the centre of the core 
at which each primary protostar forms; the shaded band 
is the transition zone defined above. Protostars are equally 
likely to form in the central zone, the transition zone, or 
the outer zone. Thus, the turbulent velocity field appears to 
erase quite effectively the initial density structure. 

THE FLOW PATTERN NEAR A CONDENSATION. Fig. |4] 

shows the positions of SPH particles in the vicinity of the 
primary protostars is Runs lb and 6d in the disc's system 
of inertia frame. The particles are colour coded according 
to their temperature. The dominant flow pattern involves 
material accreting onto the filament, then fiowing along the 
filament towards and into the protostellar disc, and finally 
ending up in the protostar. As material accretes onto the fil- 
ament it is heated to ~ 20 K (this is evidenced by the sheath 
of warm particles at the edge of the filament), and thereafter 
it is heated further by compression, reaching ~ 50 K as it 
approaches the disc. 

THE INFLUENCE OF FILAMENT DYNAMICS ON THE CON- 
DENSATION. In Run lb, the protostellar condensation is near 
the end of a filament, and so accretion onto the protostellar 
disc is quite lop-sided. In Run 6d, the protostellar conden- 
sation is close to the centre of a tumbling filament. Con- 
sequently comparable accretion fiows converge on the con- 
densation from either side. However, because the filament is 
tumbling, these flows become offset from one another (see 
Fig. [T]and the {x,y) projection in Fig. 2}. This is why the 
protostellar disc acquires a large amount of angular momen- 
tum, and is relatively extended. 



5 DISC PROPERTIES 

DISC EVOLUTION. Fig. [S] shows false-colour density images 
on the principal cartesian planes through the protostellar 
discs forming in Runs lb and 6d. We see that the disc in 
Run lb is much more compact and dense than that formed 
in Run 6d. Fig. [6] shows the protostellar mass as a function 
of time, for all the turbulent simulations. We see that in 
Run lb the protostar forms late (tg = 150 kyr), but then 
grows quite fast (M* ~ lO~^M0yr~^); this is typical of 
the majority of runs, in that the filament is not tumbling 
significantly and therefore material fiows rapidly along the 
filament and either directly into the protostar or onto a com- 
pact protostellar disc. Conversely, in Run 6d the protostar 
forms early {t^ = 112 kyr), but then grows more slowly 
{Mi, ~ 5 X lO~^M0yr~^); this is typical of a minority of 
runs, in that the filament is tumbling and therefore material 
flows initially into an extended protstellar discs and then 
spirals slowly onto the protostar. 

TERMINATION OF THE SIMULATIONS. In W09 we followed 

all simulations until 28% of the core mass was in the primary 
protostar or the protostellar disc. However, here the pri- 
mary protostar - once formed - tends to grow more quickly, 
thereby slowing down the simulation dramatically; basically 
most of the processing power is being used following the mo- 
tion of the dense material in the protostar. As a consequence 
we have had to terminate the simulations at arbitrary times, 
^END, when the evolution has become intolerably slow (see 
Table 1). 

DISC DENSITY PROFILES. Figs. [T^jb show density pro- 
files, p{lj, 2 = 0), on the equatorial planes of the protostellar 
discs produced in Runs lb and 6d. The disc in Run lb has 
p{uj,z = 0) oc Ld^^'"^, and extends out to ~ 100 AU. The 
disc in Run 6d has p{u>,z = 0) oc , inside ~ 100 AU, 
and then the profile steepens to p{uj,z = 0) cx uj~^'^ be- 
tween ~ 100 AU and ~ 1000 AU; it is more extended, and 
less dense, than the disc in Run lb. 

DISC TEMPERATURE PROFILES. Figs.[7l;,d show temper- 
ature profiles, T{w,z = 0), on the equatorial planes of the 
protostellar discs produced in Runs lb and 6d. They both 
approximate to T cx (i.e. marginally steeper than for 

the discs formed from rigidly rotating cores in W09). The 
disc in Run lb is somewhat hotter than that in Run 6d, 
because it is denser and more compact, and therefore cools 
less efficiently. 

DISC STABILITY. Fig. |S] shows the Toomre parameter, 
Q^j, , for the relatively cool and extended disc formed in Run 
6d. We see that there is a small unstable region in one of 
the spiral arms, and therefore it is possible that a fragment 
would condense out here, if the simulation could be followed 
further. All other discs seem to be gravitationally stable at 

THE ORIGINS OF PROTOSTELLAR ANGULAR MOMENTUM. 

Protostellar discs tend to be oriented with their rotation 
axes perpendicular to the filament in which they condense 
(see Fig. |4]). This implies that the net angular momen- 
tum of a filament derives from the fact that the filament 
is tumbling about one of its short axes, rather than spin- 
ning about its long axis (as, for example, in the model of 
Baneriee er aD'200( i). This a grees broadly with the finding of 
Anathpind ika & W hitworthI (2008). who analysed the pro- 
jected angles between filaments and the outfiows from pro- 
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Figure 4. SPH particles in the vicinity of the primary protostar in Run lb at time t = 168 kyr, and in Run 6d at time t = 138 kyr. 
Particles are colour coded according as their temperature is T < 20K (grey), 20K ^ T ^ 50K (orange), or T > 50K (red). The 
co-ordinate system here is centered on the protostar and rotated according to the principal axes of the disc's system of inertia. The 
circle and the two straight lines indicate the projections in which the disc is seen, respectively, face on or edge-on; the disc itself is too 
small to be resolved on this plot. The arrows give the projected direction of the disc angular momentum. The arrows are normalised to 
j = 1.03AU^s~^ for Run lb, and to j = 0.19AU^s~l for Run 6d, respectively. 




Figure 8. False-colour image of the Toomre parameter, Q^, for 
the disc formed in Run 6d, at time t = 144.8 kyr. The centre of 
co-ordinates has been shifted to the centre of mass of the disc. 



tostars embedded within the filaments. On the assumptions 
(a) that the outflow direction corresponds to the angular mo- 
mentum of the underlying accretion disc, and (b) that their 
sample was randomly oriented and not subject to any selec- 
tion efl'ects, they inferred that most outflows are at a large 
angle {9 > 45 ) to the fllament in which the driving proto- 
star is embedded. However, w e should emphasise that the 
observed filaments analysed by I Anathpindika fc WhitworthI 
are larger than the prestellar cores they contain, whereas 
the filaments formed here are within prestellar cores. 

DISC PRECESSION. From Fig. |4]we see that the angular 
momenta of the protostellar discs are not in general aligned 
perpendicular to the instantaneous plane of the disc. This in- 
dicates that the discs are precessing, which is to be expected, 
because, in a turbulent environment, the disc is assembled 
by irregular, asymmetric and lumpy inflow, and the direc- 
tion of its net angular momentum is repeatedly changing. 



6 CONCLUSIONS 

We have simulated star formation in prestellar cores with 
transonic turbulence. The cores are all modeled as super- 
critical Bonnor-Ebert spheres (^3 — 6.9), with the den- 
sity increased by 10%. All the cores have the same mass. 
Mo = 6.IM0, and the same radius, = 17, OOOAU, but 
different realisations of the initial turbulent velocity field. 
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Figure 5. False colour density slices on the principal cartesian planes through the protostellar discs formed in Runs lb and 6d, 
with velocity vectors superimposed. The three lefthand columns show frames from Run lb; successive rows correspond to the times 
t = 148, 150, 152, 154, 162 and 168 kyr. The three righthand columns show frames from Run 6d; successive rows correspond to the times 
t = 112, 113, 119, 129, 137 and 145 kyr. The density scale is shown at the top of the page, and a lkms~^ velocity vector is given on the 
a; = slices. All frames are (1000 AU)^, and the centre of co-ordinates has been shifted to the centre of mass of the disc. 
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Figure 7. Density and temperature profiles on the equatorial planes of the protostellar discs formed in Runs lb and 6d. (a) p{u!, z = 0) 
for Run lb. (b) p(t<j, 2 = 0) for Run 6d. (c) T{ui, z = 0) for Run lb. (d) T{uj, 2 = 0) for Run 6d. The centre of co-ordinates has been 
shifted to the centre of mass of the disc. 



The simulated realisations have been selected from a large 
ensemble of four hundred different realisations, so as to span 
the same range of specific angular momenta, j, as the simula- 
tions performed in W09. In the rigidly rotating cores treated 
in W09, a collapse wave propagates from the outside in, and 
the density in the centre of the core increases monotonically 
until this collapse wave arrives after about one freefall time, 
tpp. At that stage a protostar forms, surrounded by a pro- 
tostellar disc. The extent of the disc depends on j, and only 
cores with large j produce extended protostellar discs which 
are prone to fragment. In contrast, the key features of the 
turbulent cores simulated here are the following. 

• The outcome is not strongly dependent on j. Rather the 
outcome reflects the stochastic nature of the initial turbulent 
velocity field. The gas at the centre of the core is perturbed 
from the outset, and turbulent fiows tend to sweep up the 
gas of the core into a filament, whose extent is comparable 
with the diameter of the core (~ 30, 000 AU). 

• Protostars condense out of these filaments. If two proto- 
fragments form in a filament at widely separated locations, 
they may survive to form a wide binary system, but usually 
proto-fragments in the same filament tend to merge before 
they can condense out. Proto-fragments can form at loca- 
tions significantly displaced from the centre of the core, and 
they form with peculiar velocities of order 0.1 kms~^. 

• Proto-fragments condense to form a protostar and then 
an attendant protostellar disc. However, these protostellar 
discs tend to be quite compact, with a shallow density pro- 
file, and stable against gravitational fragmentation or the 
formation of spiral waves. Part of the reason why they are 



stable appears to derive from the extra dynamical heating 
they experience as a result of lumpy, irregular accretion from 
the turbulent core envelope. 

• The intrinsic angular momentum of a proto-stellar disc, 
and hence also its rotation axis, tend to be perpendicular to 
the long axis of the birth filament. This suggests that the 
angular momentum of a proto- fragment derives from the fact 
that the filaments are tumbling, rather than spinning about 
their long-axis. 

• A key parameter infiuencing the statistics of turbu- 
lence in a core is the minimum wavenumber, fc^jj., . Only 
if fcj^jtj = 1 (i.e. turbulence is injected on the scale of the 
diameter of the core) do we reproduce the range of spe- 
cific angular momenta observed in nature. Thus it is the 
largest wavelengths in the turbulent spectrum which are 
critical. This accords with the result reported recently by 
iBatd (|2009bl ). who shows that the statistics of stars formed 
in clouds with oc fc^* and Pk oc are indistinguishable. 
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